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Abstract 

The purpose of this paper is to explain the discrepancies existing in 
the literature relative to e^e~ pair production in peripheral heavy ion col- 
lisions at ultra-relativistic energies. A controversial issue is the possible 
cancellation of Coulomb corrections to the Born term in the pair pro- 
duction cross-section. Such a cancellation has been observed in a recent 
approach based on finding retarded solutions of the Dirac equation, but 
does not seem to hold in a perturbative approach. We show in this pa- 
per that the two approaches are in fact calculating different observables: 
the perturbative approach gives the exclusive cross-section of single pair 
production, while the other method gives the inclusive cross-section. 

We have also performed a thorough study of the electron propagator 
in the non-static background field of the two nuclei, the conclusion of 
which is that the retarded propagator is in the ultra-relativistic limit a 
much simpler object than the Feynman propagator, and can be calculated 
exactly. 

PACS codes: 25.75.-q, 11.80.-m, 34.90. -hq, 11.55.-m BNL-NT-01/1 
Keywords: Ultra-relativistic nuclei collisions, Electromagnetic pair production, 
Eikonal approximation. Coulomb corrections. 



1 Introduction 

In the past three years, the problem of e^e~ pair production induced by the 
collision of two nuclei has attracted a lot of interest, due partly to a series of 
papers showing that the Dirac equation of an electron can be solved exactly 
in the electromagnetic background field created by the two nuclei, in the limit 
where the two nuclei are ultra-relativistic j^, |2| || ^ . This solution was then used 
to derive an expression for the pair production cross-section. An unexpected 
consequence of this formula was that it lead to a pair production cross-section 
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equal to its Born value (given by the diagram in figure In other words, this 
result seems to indicate that all the corrections due to the multiple exchange of 
photons (called "Coulomb corrections" in the following) cancel exactly in the 
total cross-section. 

< i < 
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> ? > 

Figure 1: Born contribution to the pair production amplitude. 

This observation then prompted another series of papers on the same topic, 
basically saying that this property cannot be true. One of them |^ re-derived 
with new methods the classical result of Davies, Bethe and Maximon |^, || 
about the fragmentation of a photon into a pair, induced by a electromagnetic 
background field. Then, this result was applied to pair production in nuclear 
collisions ^, In this perturbative calculation, the Coulomb corrections do 
not cancel and the resummed cross-section is about 25% smaller than the Born 
cross-section at RHIC energy. Others claimed the discrepancy to be due to a 
breakdown of crossing symmetry Another paper attributed the con- 
tradiction to a mishandling of ill-defined integrals, and claimed that the result 
of Ivanov et al. can be recovered from the solution of the Dirac equation if a 
proper regularization is used. But those papers apparently focused on techni- 
cal details, without first asking if an agreement between ||, |[ Q and |9[ ^ 
should be expected at all. In other words, are those two approaches calculating 
the same physical quantity ? Given the fact that one obtains a retarded ampli- 
tude by solving the Dirac equation, while the perturbative approach calculates 
a time-ordered amplitude, such an assumption is highly non-trivial, and in fact 
very unlikely. 

In this paper, we re-examine from first principles the problem of pair pro- 
duction by two colliding nuclei. Our main focus is on showing that the two 
approaches used to study pair production in the peripheral collision of two 
ultra-relativistic nuclei are in fact calculating different physical quantities. In 
particular, we give the relationship between several physical quantities and the 
appropriate correlators. The exclusive cross-section of single pair production is 
related to a time-ordered (or Feynman) correlator, and is the object usually cal- 
culated in perturbation theory. On the other hand, the inclusive cross-section 
of pair production can be related to retarded amplitudes, which are the ones 
accessible by solving the Dirac equation with retarded boundary conditions. 
Therefore, the main discrepancy stems from the fact that the two approaches 
are calculating different quantities, despite giving them the same name; and the 
two approaches do not have to give the same answer. 

We then pursue this analysis by studying the retarded and Feynman prop- 
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agators in the background field of two ultra-relativistic nuclei, and show that 
these two types of propagators have dramatically different perturbative expan- 
sions: the retarded propagator can be obtained in closed form in the ultra- 
relativistic limit, as opposed to the Feynman propagator. 

We also resolve a puzzle related to the unitarity of the pair production 
probabilities calculated in perturbation theory, which seems to give probabilities 
larger than 1 |l^, |l^, We show that this calculation must include as 
a corrective factor the vacuum-to-vacuum transition probability, and that this 
factor precisely restores unitarity. 

The structure of this paper is as follows: after setting up some notations 
and details about the model in section ||, we recall in section ^ some relations 
between observable quantities and various types of correlators. We also point 
out some uncommon aspects of the Feynman rules that are due to the fact that 
the background field is non-static. 

We study in section ^ the propagators in the background field of only one 
nucleus. This case can be solved exactly in the ultra-relativistic limit, and serves 
as a starting point for the more complicated case of two nuclei. 

In section ^, we begin with a formal solution to the problem of two nuclei, 
which enables one to write both the retarded and Feynman propagators in terms 
of the scattering matrices associated to the individual nuclei. Then, we explain 
why the ultra-relativistic limit simplifies this formula in the case of the retarded 
propagator, but does not help for the Feynman one. 

In section ^ we focus on the problem of the unitarity of the perturbative 
predictions, and show how it is resolved. We also show that the multiplicity 
distribution is not exactly Poissonian. 

Section |^ discusses two strategies one can use to get an approximate answer 
for the Feynman propagator, and hence for the pair production amplitude. 

Finally, section p| is devoted to concluding remarks. 



2 Ultra-relativistic limit. Notations 

In this paper, we consider two nuclei of electric charge Z\e and Z^e respec- 
tively, colliding with an impact parameter h. We study the problem in the 
frame where the two nuclei have opposite velocities. For the sake of definite- 
ness, we decide that the nucleus Z\ is moving in the positive z direction, while 
the other nucleus is moving in the negative z direction. We denote by A^i^.x) 
and ^2 (t, x) the electromagnetic vector potentials created by those nuclei. Be- 
cause the superposition theorem holds in QED|^, the total vector potential is 
A^'i^t, x) = A^{t, x) + At^{t, x). 

In the ultra-relativistic limit, we can assume that the two nuclei do not recoil, 
and simply treat them as a classical background field. The trajectories of the 
two nuclei are respectively z = ztt {c — 1), x± — ±6/2, and there is a gauge in 

^This property is not true in QCD. The problem of qq pair production is therefore compU- 
cated by the necessity of first finding the color vector potential due to the two nuclei. 
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which A^l and A2 have the following form |17, 19 : 



A'^{t,x) = ZieS{v+ ■ x)v'^ In 
At^{t,x) = Z2e(5(u_ -x)!;^ In 




(1) 



where v'^ = (1, 0, 0, ±l)/\/2. This choice of gauge makes the following fact 
obvious: the electromagnetic field of an ultra-relativistic nucleus is confined in 
its transverse plane due to Lorentz contraction. In the ultra-relativistic ap- 
proximation, we neglect any modification of the trajectories of the two nuclei. 
Therefore, one can see this problem as a field theory for fermions in a classical 
electromagnetic background field. An essential property of the field generated 
by two colliding nuclei is that there is no frame in which it is time-independent 
(static). This property has important consequences on the field theory describ- 
ing fermions in this background; one of them is that pairs can be created^. 

However, one should realize that the limit of infinite momentum for the nuclei 
is only an approximation of the real physical situation. Indeed, the potentials of 
Eq. (|l|) cannot be expected to describe the motion of a lepton comoving with one 
of the nuclei. This restriction also appears in the work of |2^ . In particular, 
shows that the gauge transformation that leads to Eq. (|l^ is well-defined 
only if the lepton is not comoving with one of the nuclei. There is another 
way to see the problem with Eq. (|l|) : this background potential is short-ranged 
in the z and t directions, while the potential before the gauge transformation 
was long-ranged. The absence of interaction at asymptotic times allows the 
construction of plane- wave asymptotic states, and makes the usual formulation 
of reduction theory applicable. But it is also clear that those interactions cannot 
be removed by a gauge transformation for a comoving lepton, the states of which 
are described by distorted waves. In conclusion, those potentials are valid only 
if the typical interaction time between the lepton and the nuclei goes like I/7, 
i.e. for leptons produced in the mid-rapidity region. 

In addition, the description of the background field by a delta function is 
probably not accurate at large transverse distances if the Lorentz factor 7 is 
finite. Indeed, at large distances, the lines of force of the electric field should 
have some curvature and depart from the transverse plane. Therefore, Eq. 
is expected to break down also at large impact parameter b or if the transverse 
separation between the lepton and a nucleus is large. This problem is expected 
to show up via infrared divergences in the transverse momentum integrals. In 
[||, it was argued that those divergences are regularized by an effective cutoff 
of order cu/^. Although this argument cannot tell precisely what this cutoff 
should be, it is obvious from the derivation of the ultra-relativistic limit of the 
potentials that it is completely determined by the kinematics and is therefore 

■^This is to be contrasted with the case of a single nucleus. In the frame of that nucleus, 
the background field is static, and it is well know that a single nucleus moving on a straight 
line does not produce pairs (except via the non-perturbative Schwinger mechanism of vacuum 
instability |20|, which we ignore here). 
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independent of ^1^20; (the atomic numbers and coupling constants appear only 
in the numerator of the Coulomb potential). Therefore, the uncertainty in the 
cutoff can affect only the overall normalization of integrated cross-sections, but 
not its dependence on atomic numbers. 

In view of those caveats, one should take with a grain of salt the results 
obtained in this model for completely integrated (over the pair phase-space and 
over impact parameter) cross-sections, because there are regions in such an 
integral where the above approximation may not be appropriate. 

One should realize that due to the large electric charge Ze of a nucleus like 
gold, the addition of an extra photon connecting the electron line to a nucleus 
brings a factor Za (a = e^/47r), which may be of order 1. Therefore, we want 
to include as much as possible of these "perturbative" corrections. On the 
contrary, photons coupling an electron line to itself (or to another electron line) 
yield only a factor of a, which is indeed a small correction. Therefore, we also 
neglect the interactions of the electrons with dynamical photons, and keep only 
the classical photon background. The Lagrangian density for the electron field 
'^(x) is therefore 

£ = '^{x) {i0^ ~ e4(x) - m)4>{x) (2) 

in this model. 

In the following, we make extensive use of the light-cone coordinates. For 
any 4- vector x^ , we define: 



„± — 



x^ ± x^ 



and denote by x± the transverse part of the 3-vector x. With these notations, 
the invariant norm of x^ is x^ = 2x'^x~ — x±^, and the scalar product of fc^ 
and x^ is k ■ X = k^x" + k~x^ — k± ■ x±. The invariant measure d'^x becomes 
d'^x = dx~^dx^d^x±. Note also that x^ ^ x ■ Vzfz with the v± defined above. 



3 Reduction formulae and Feynman rules 

In this section, we relate observable quantities to correlation functions of the 
fermionic field operator. In fact, these considerations do not depend on the 
nature of the background field, and rely only on the Lagrangian of Eq. (||) . In 
particular, the formulae of this section are independent of the ultra-relativistic 
approximation. 

3.1 Amplitude to produce one pair 
3.1.1 Reduction formula 

The amplitude to produce one e^e~ pair is 

Mi(p,g) = (e+(p)e-(q)out|Oin) = (Oout|dout(p)&out(q)|0in) , (4) 
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where rfout(p) (resp. &out(p) ) is the operator that creates a positron (resp. an 
electron) of S-momentum p in the final state. At this stage, it is very important 
to carefully distinguish the in- and out- states and operatorg|^. 

Making use of the following relations between the annihilation/creation op- 
erators^ and the field itself (see page 61): 

boutiq) = J d^a; ?l(q)7>out(i,a7)e'«-^ , 

doutip) = J d'x ^,^,{t, x)^%ip)^P-- , (5) 

where po and go are, respectively, [p^ + w?) and [q^ + m?) (the time = t 
that shows up in these formulae is irrelevant, and disappears in the course of 
the calculation), one can show by standard manipulations ([^, pages 205-207) 
that the pair production amplitude defined above can be related to a 2-point 
time-ordered correlator by the following reduction formula: 

■ n 2 



(0outMout(p)&out(9)|0in) = 



d'^x d'^y 



xe''>--u{q){i 0., -m) (Oout|T^(y)V'(x)|Oin) (* 0y +mHp)e'P-y , (6) 

where the arrows indicate on which side the derivatives act (preventing their 
action on the exponentials), and where Z2 is the wave function renormalization 
factor for an electron^. The important point to note here is that the average 
value of the time-ordered product is taken between the in- and out- vacua, which 
are different states. To expand a little on this, let us add that we necessarily 
have I (Oout|0in) P < 1 in background fields that can produce pairs, because 
of unitarity (see section || for the role played by (Oout|0in) in issues related to 
unitarity) . 

3.1.2 Perturbative expansion of (Oout|T-!/)(j/)'(/;(a;)|Oin) 

We need to calculate the correlator (Oout|T 7/;(y)V'(a;)|0in). From that, Eq. (||) 
tells us how to obtain the pair production amplitude: amputate the external 
legs of the correlator, take its Fourier transform, and insert the result between 
the appropriate spinors. 

The Feynman rules to calculate perturbatively this correlator can be ob- 
tained by switching to the interaction picture. The Heisenberg field can be 
expressed in terms of the field in the interaction representation, via the follow- 
ing relation: 

il;{x) = U{t,,x°)'iJj,{x)U{x°,t,) , (7) 

•^For instance, feout does not annihilate the in- vacuum: wout |Oin> ^0, while feoutlOout) =0. 
It is precisely this fact that makes the pair production amplitude non-zero in a non-static 
background field. 

*For the sake of brevity, we do not write explicitly the spin indices on spinors and cre- 
ation/annihilation operators. In formulae for cross-sections, it is implicitly assumed that we 
sum over the spin of the final state particles. 

^With the Lagrangian of Eq. (H), this factor is equal to 1. 
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where is the time at which the Heisenberg and interaction pictures coincide 
(ultimately, we will take to —oo), and where the evolution operator U is 
related to the interaction part C-mt = —etp{x)4-{x)ip{x) of the Lagrangian by 

U{t2,ti) = Pi2expi f d'^x Cintii^,{x)) , (8) 

where P12 is an ordering operator along the path connecting ti to ^2 (ordinary 
time-ordering if ti < ^2, and reverse time-ordering if ti > 12)- Using this 
transformation, and taking —00 (in this limit, tp^ — > '>jJin), we have 



(Oout|T?/'(y)V'(a;)|Oi„) = ^Oout Ui-oo,+oo)Tip,^{y)iPi4x) 

r+oo 

Xexpi/ d'^x£int{lpin{x}) OinJ ■ (9) 



Noticing that (Oout|^^(— 00, +00) = (Oi„| , we have simply: 
(0out|T7/i(y)V'(a;)|0in) = (Oin T7/;i„(y)?/;in(a;)expi / d'^x£int(V'in(a;)) 



Oi, 
(10) 

Everything being now expressed in terms of in-fields and in-states, the right 
hand side of the last equation can be evaluated by the standard perturbative 
expansion of the exponential. 

This perturbative expansion has, however, a peculiarity due to the fact that 
the background field is non-static: the vacuum-vacuum diagrams (i.e. diagrams 
without any external legs) do not cancel, and their sum is not a phase. On the 
contrary, in a conventional field theory, (Oout|0in) is just an irrelevant phase, and 
one takes advantage of this fact to divide the r.h.s. of the previous equation 
by (Oin|[^(+cx3, — cx3)|0in). Then, one can show that this denominator cancels 
( |2^ , pages 266-267) order by order the vacuum- vacuum diagram^ This trick 
cannot be used here due to the fact that (Oout|0in) is not a phase in the present 
problem ||2^ . In fact, this is deeply rooted in the property that the background 
can produce particles. 

We can keep this complication aside for a while by just writing: 



(0out|TV'(2/)V'(2;)|0i„) - (Oin|C/(+^,-C»)|Oi, 



Oi, 



T'0in(2/)V'in(a;) exp i f^^ d'^xCint{'4'in{x)) 



Oi, 



(0i„|?7(-fCX),-CX))|0in) ' ^^^^ 

so that the fraction on the right hand side has a perturbative expansion where 
the vacuum-vacuum diagrams do cancel. In fact, this fraction is nothing but 
the Feynman propagator 

(Oout|TV;(y)^(x)|Oin) 



''This step is not mandatory though. One can live with the vacuum-vacuum diagrams and 
notice that they add up to a pure phase, so that they always drop out in the calculation of 
cross sections, even if they are present in the amplitude. 
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of an electron in the electromagnetic background field. Its perturbative expan- 
sion is the usual one: it is obtained by inserting the external potential —ie4^{x) 
on chains of free Feynman propagators G°^{x,y) = (Oin|T?/;jj^(y)^/;in(a;)|Oin)- In 
particular, the perturbative expansion for the full propagator can be generated 
by the following Lippmann-Schwinger equation: 

G{x, y) = G°{x, y)-ie [ d'^z G°(x, z)4{z)G{z, y) . (13) 



Note that this equation is equally valid for the Feynman and for the retarded 
propagator (provided the free propagator G^ is chosen accordingly). This re- 
mark will become important later when we also need to study the retarded 
propagator. 

It turns out to be more convenient to work in Fourier space. If we define[| 

" / (0J (0^^""''^""^('^'^) ' 
then the Lippmann-Schwinger equation becomes 

G(q,p) = {2i,)H{p - q)G"{p) - leG'iq) f ^A{k)G{q + k,p), (15) 



{27:y 

with G'^{p) = ijij^^m) (and an unspecified it prescription depending on whether 
we are studying the retarded or the Feynman propagator) and 

A^{-k) = [ d^x e'^-''A^'{x) . (16) 



If one introduces the interacting part Tp of the Feynman propagator by the 
relation 

GF{q,p) = {2t:)H{p - q)G'i.{p) + G'i.{q)TF{q,p)G'^^{p) , (17) 

then Eqs. (^ and ( pi] ) lead to the following expression for the probability to 
produce exactly one pair in a collision at impact parameter b: 

Note that this formula has been derived independently of the details of the 
background field, and is therefore completely general. 

3.2 Average number of produced pairs 



^We do not use a distinct symbol for Fourier transforms, as the context always enables to 
tell the difference. 
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3.2.1 Expression as a correlator 



The approach based on solving the Dirac equation was motivated by the papers 
[ p5| p6| , [27|] , where a formula giving the average number of pairs produced per 
collision is derived in terms of retarded amplitudes only. We present here a 
justification of this formula in the ficld-thcoretical framework we are following 
in this paper. 

Let us start from the expression of the average number of pairs n as the sum 



+00 



+00 



n = ^ nPn + l)Pn+i , 



(19) 



ri=0 



where we denote by P„ the probability to produce exactly n pairs in a collision 
at impact parameter b. Making explicit what this probability is, we find 



n - 



+00 
x|(0o 



1 



n+1 

n 



dh 



2uJa 



dout(Pl)&out(qi) ■ • ■ rfout(Pn+l)&out(9n+l) 



Oi, 



(20) 



where we denote lo^ = ^/(p^ + m}). Singling out one of the electrons (say the 
one with momentum which we call simply q) and expanding the squared 

modulus, we can write: 



d3q 



Oi, 



^Lt(g) 



V f TT ^ TT ^ 

n! (n+1)! ] ^\ (27r)32a;p, ^\ (27r)32c^q^. 

n— i— 1 ^ J — 1 ^-^ 

4ut(Pl)^Lt(9l) • • • 4ut(Pn)&Lt(9n)4ut(Pn+l) 



Oout 

x6out(g) 



rfout(Pl)&out(gi) ■ • • rfout(P«)&out(g«)rfout(Pn+l^ 



(21) 



Noticing now that the three intermediate lines are the identity operator on the 
subspace of states with electric charge +e, we find: 



(27r)32wq 



Oi, 



&Lt(g)^out(g) 



Oi, 



(22) 



This formula simply tells that in order to count the number of pairs produced if 
the initial state is the vacuum, it is sufficient to count the number of electrons 
in the final state. Had we decided to single out a positron annihilation operator 
in Eq. ([2l]), we would have obtained instead: 



d^p 



(27r)32wp 



Oi, 



<it(P)'^out(p) 



Oi, 



(23) 
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One can note that if the background potential has time-reversal symmetry, we 
have: 



(27r)32wq 



bl{q)bin{q) 



(27r)32cjp 



d\^(p)din(jp) 



(24) 

There are also reduction formulae for those correlators, that give for instance 



d^q 



(27r)32w, 



d xd y 



xe''i--u{q)[i 0^ -m) {Q,^\'^{y)iP{x)\K-) {i 0y -m)u{q)e-^'i-y . (25) 

Therefore, in order to calculate the average number of pairs, we need now the 
ordinary product of two fields, averaged with the initial vacuum |Oin) on both 
sides. 

3.2.2 Perturbative expansion of (Oin|'0(j/)V'(2;)|Oin) 

In order to switch to the interaction representation, it is easier to start with 
(Oin|T'0(y)?/'(x)|Oin) for which one can write directly 



(0i„|TV'(2/)V(a;)|0i„) 



[/(-cx), +cx))T-0ij^(?/)-0in(a;) 

r4. 



X exp I 



d xCintiAnix)) 



Oi, 



(26) 



This time, one cannot get rid of the C/(— oo, +oo) in the right hand side, but there 
is a standard trick to incorporate it in the perturbative expansion j28| ^ ^ . 
For that, one has to introduce a contour C going from — oo to -\-oo just above 
the real-time axis, and then back from +oo to — oo below the real axis: 



One can then check that 



(Oi„|Ti^(y)^(a:)|Oin) = (Oi„ ¥^,^{y)t.n{x) expi / d xCint{i^in{x)) 



(27) 



In the previous formula, P stands for an ordering of operators along the path 
C, identical to the usual time-ordering T on the upper branch of the contour 
(and with the convention that points on the lower branch have a "larger" time 
than points on the upper branch) . We therefore have formally similar Feynman 
rules for the perturbative expansion of (Oin|Ti/)(y)t/;(x)|Oin), except that the 
time integrations at each insertion of the external potential are performed on 
the contour C. 
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It is customary to write this formalism in matrix form | p8| , |3l| (in this 
section, objects denoted by a boldface letter are 2x2 matrices), by splitting 
the free propagator in four components according to where and y° lie on C 
(upper or lower branch): 

GO+(x,y) Gl_{x,y) 



with 



Gl+{x,y) = (Oi„|TV^i„(y)V'i„(a;)|Oi„) , 

G°„(x,2;) = (Oin|T^i„(y)^i„(a;)|Oi„) , 

G°_(x,y) = (0i„|^i„(2/)Vin(a:)|0in) , 

G%{x, y) = (Oi„| - ^i„(a;)^i„(y)|Oi„) , (29) 

where T is the reverse time-ordering operator. The correlator (Oin|'0(2/)'0(2;)|Oin) 
we need for the average number of produced electrons is the H — component of 
the exact matrix propagator. 

One needs also to give a matrix structure to the external potential^(a;): 

4{x) = T34{x) , (30) 

where T3 = Diag(l,— 1) is the third Pauli matrix^ Then, at each vertex, the 
rule is to integrate only over the ordinary time axis, and to multiply the ma- 
trices corresponding to the propagators and potential insertions (in the order 
they appear in the Feynman diagram). For instance, the Lippmann-Schwinger 
equation for the exact matrix propagator in Fourier space is 

G{q,p) = {2nf5{p - q)G\p) - ieG\q) J -^^j^{k)G{q + k,p). (31) 

Note that in Fourier space, the expression of the free matrix propagator is^: 

gUp)-^ J^T^- =[g°-(p)]% 

— m'^ + le 

Gl_{p) = 27r0i-poW+m)S{p^ - m^) , 

G"_+(p) = 27r0(+po)(y+ m)5(p2 - m^) . (32) 

As it stands, the expansion of Eq. (^) has a very intricate matrix struc- 
ture. However, it can be simplified by applying a "rotation" ^ on the 

*The — sign for the component of^ comes from the fact that the integral over the 

lower branch of the time contour goes from +00 back to — Qfl, 

^The expression of those propagators can be found in |p2|l (Eqs. (3.92)3. To apply them 
here, the following substitutions must be made: (7 = 0, n(po) = 0. Also, |32| is using different 
notations for the indices labelling the two branches of the contour: 1 = +, 2 = — . 
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previously defined matrices. Let us define 

G'ij{p)^U{p)G'{p)U^{-p), 
Gu{q,p) = U{q)G{q,p)u"{-p), 
4u{q, q + k) = u"-\-q)4{k)U-\q + k) 



(33) 



where U is an invertible matrix. The Lippmann-Schwinger equation hardly 
changes in the rotated formalism, 

Gu{q,p) = {2T,)H{p-q)G'lj(p)-ieGl{q) J ^4u{q,q + k)Gu{q + k,p) , 

(34) 

but there are some choices of U that simplify significantly the free matrix prop- 
agator. A convenient choice is 



which leads to the free propagator 



-1 



GVrt 

with the definitions (the index "S"' stands for "on-shell"): 



(35) 



(36) 



G%{p) = 2TT{^+m)5{p-'-m^) . 



jp — m? — ipoe 



In this transformation, the external potential becomes 

4uiq,q + k) =4ik) ( J J 



(37) 



(38) 



Therefore, we have a simplification in the propagator which has now a vanishing 
component, the non-zero components being the free retarded and advanced 
propagators and the very simple on-shell piece G^. That this transformation 
helps to rcsum the pcrturbativc expansion is readily seen by computing the 
building block that will be iterated in the expansion: 



4u{q,q + k)Glj{q + k)=4{k) 



G%{q + k) 
C^iq + k) 



+GUq + k) 



1 




(39) 

Indeed, the fact that this object is the sum of a diagonal matrix and a nilpotent 
matrix leads trivially to the following formula for a term with n insertions of^: 



Guiqo) 4u{qo,qi)G°i;{qi) ■■■ 4u{qn-i,qn)G%{qn) 
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GO (9o)4(gi - go)Gj\(qi) • • 4{qn - qn-l)G\{qn) 



- J2 G'R{qo)4{qi - go) • • • G%iq,^Miq^ - q^-l)GUq^ 



x4(g-*+l - g^)G5l(q^+l)4(q,+2 



q^+l) 



■G\{qn) 




1 



(40) 



At this stage, it is straightforward to sum over the number n of insertions of the 
external potential, to obtain the exact propagator in the new basis (set q = qq 
and p = Qn in the previous formula) : 



Guiq,p) 



GA{q,p) 
GR{q,p) Gs{q,p) 



(41) 



where GR^A{q,p) are the exact retarded and advanced propagators, and where 



Gs{q,p) 



-2Ti5{k^ - m^) 



(2^)4 

xGfl(<z, k)Gl-\k){H + m)G\-\k)GA{k,p) 



(42) 



At this stage, we have an exact expression for the Fourier transform G^ {q,p) 

of (0in|''/'(2/)V-'(2;)|0i„), which reads 



G+-{q,p) = - GAiq.p) - GR{q,p) + Gs{q,p) 



(43) 



Using Eq. (|4^), we see that this correlator can be written in closed form in 
terms of Gr and G/tP]. According to Eq. (^), n is obtained by amputating 

the external legs of G-| (q, g). To that effect, it is convenient to introduce 

the interacting parts TrIq^p) and TAiq^p) of the exact retarded and advanced 
propagators, defined by 



GR{q,p) = {2^)H{p-q)Gl{p)+Gl{q)rR{q,p)Gl{p) , 
GA{q.p) = {2T,f6{p-q)G\{p) + G\{q)TA{q.p)G\{p) , 



(44) 



which leads to 



-u{q) TA{q,q) - TR{q,q) 



(27r)22cj5 

^27r<5(fc2 _ vn?)TR{q, 



+ m)rA{k,q) u{q) . (45) 



^'^It is now clear why Eq. ([jl]) had a very compHcated matrix structure: trying to solve 
directly this matrix equation amounts to write the exact retarded and advanced propagators 
in terms of bare time-ordered propagators, which is not easy in practice. This explains why 
the rotation leading to the retarded/advanced basis simplified a lot this calculation. 
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By using the Lippmann-Schwinger equations for the retarded and advanced 
propagators, one can check thatP] 

xTR{q,k){H + m)rA{k,p) . (46) 
From there, we can further simphfy n, and obtain: 

<pQ r d^k 



Xu{q)TR{q, k)m + m)TA{k, q)u{q) . (47) 

Changing k ^ —k in the above formula, and integrating over fco thanks to 
the 6{k^ — rn?) permits to simphfy even more the expression of n. Noticing that 

TA{~k,q) = -[TR(q,-k)\\ (48) 

- m ^^v{k)v{k) , (49) 

spin 

we finaUy obtain the following compact expression for the average number of 
pairs produced per collision at impact parameter b: 

"^'^ ^ "^'^ u{q)TR{q,^p)v{p)\ (50) 



(27r)32wq J {2ttY2uj. 

Therefore, we have obtained for the average number of produced pairs a formula 
similar to Eq. (|l8|), with however two major differences: the factor | (Oout|0in) P 
is not present here, and the formula for n involves the retarded propagator 
instead of the Feynman propagator. 

3.3 Discussion of earlier literature 

Equations (|l^) and (|5^) show clearly the main difference between the approach 
based on the Dirac equation and the perturbative approach concerning the "pair 
production cross-section". Indeed, what is called "cross-section" in ^ || ^: 

CTT = J (fbn, (51) 

is the inclusive cross-section^ of pair production (obtained by counting all the 
pairs produced); while according to the nature of the diagrams considered in 
[|[ |o|, the object they called "cross-section" is in fact|^ 

a,= f (fb Pi , (52) 



^^This relation can be seen as a form of the optical theorem. 

^^Integrating over the entire range of impact parameters might include contributions from 
a region where our model is not expected to be valid. Cutoffs are implicitly understood in 
this integral, in order to make it finite. 

"'^''Up to the factor | (Oout|0in) P, which they seem to have overlooked. 
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i.e. the exclusive cross-section to produce exactly one pair (measured by count- 
ing only those events that produce exactly one pair). It is now obvious that 
one should not expect an agreement between the two approaches, since they are 
calculating different observables. In fact, it is obvious that (j\ < (Jt if there are 
collisions that can produce more than one pair. We are now in a position to 
discuss in detail the existing literature concerning the discrepancy between the 
two approaches. 

Papers @, § | 

It seems that those papers overlooked the results of ^ The latter 

papers, as well as the above derivation, only connect squares of retarded am- 
plitudes to average numbers of particles (or other moments of the multiplicity 
distribution) , but cannot give the pair production amplitude in terms of retarded 
solutions of the Dirac equation. The present derivation shows that [|l], ^, ^ are 
not correct when they call u{q)Tji{q, —p)v{p) the pair production amplitude. 
The only use of this object is via its square in Eq. (|5^), where it leads to n. 



Papers g P 

Eichmann and collaborators suggested in that the discrepancy between the 
two approaches was a consequence of the fact that crossing symmetry is not 
valid in the ultra-relativistic limit. Here is how their argument goes: in Q, they 
claim that solving the Dirac equation in the background field of the two nuclei 
leads to the exact scattering amplitude of a lepton by the nuclei. From there, 
one would have to use crossing symmetry (i.e. change p —p for the incoming 
electron) in order to obtain the pair production amplitude, a procedure which 
they later claimed to be incorrect because crossing symmetry does not work 
when the nuclei are ultra-relativistic ||Tl|| . 

However, the premises of this explanation are incorrect, because solving 
the Dirac equation with retarded boundary conditions does not give the exact 
scattering amplitude in a relativistic theory if the background potential is time- 
dependent. Indeed, in a relativistic field theory the free retarded and Feynman 
propagators differ by their ie prescription: 



Gr{p) = I— J— r , G'p{p) = I— . (53) 



The two propagators will lead to equivalent results only if the background po- 
tential cannot change the sign of the energy po (the energy of the incoming 
lepton is of course positive), i.e. only if the background potential is static. 
But if the background potential has some time dependence, it can change the 
sign of the energy carried by the propagator, and the retarded and Feynman 
propagators will lead to different results. Basic requirements of any field the- 
ory (Lorentz covariance, unitarity and causality) imply (see fs^ , pages 197-220) 
that scattering amplitudes are given by the time-ordered propagator, and not 
by the retarded one. Physically, the Feynman propagator takes into account 
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the fact that an external field may create or annihilate pairs of leptons |Q , an 
effect which affects lepton scattering. This effect is not taken into account if 
one uses the retarded propagator^ Of course, if the external field is static, or 
if the problem is treated non-relativistically^, pair creation/annihilation is not 
possible, and both retarded and Feynman propagators give the same scattering 
amplitude. 

The authors of Q also proposed an alternate proof in section III for their 
expression of the scattering amplitude, this time by a direct calculation of Feyn- 
man diagrams. However, the diagrams displayed in the figure 3 of reference 
do not vanish in general if calculated with Feynman propagators. Although it 
happens that the diagram with the configuration ABA (see reference Q for the 
notations) of external potentials vanishes, the diagram ABBA for instance does 
not. 

In conclusion, solving exactly the Dirac equation does not give the exact 
lepton scattering amplitude, which makes the discussion of crossing symme- 
try irrelevant for the present problem. Indeed, solving the Dirac equation (or, 
equivalently, calculating the retarded propagator), gives directly n thanks to 
Eq. (pO|). Of course, since one needs Tfi{q,—p) and not TR(g,p), one should 
be careful not to assume that the incoming energy is positive when calculating 
Tr. In the next two sections of the present paper, we discuss in more detail the 
differences between the retarded and Feynman propagators in the background 
field of the two nuclei. In particular, a closed form expression for Gji{q,p) is 
derived (which is valid for energies of any sign, and therefore can be used in 
Eq. (pO|)). On the other hand, we show that it is not possible to obtain such a 
simple expression for the Feynman propagator. 



Paper 

In their paper , Lee and Milstein proposed a different explanation for the dif- 
ference between the two approaches. Starting from a formula equivalent to our 
Eq. (^l]), and manipulating possibly ill-defined integrals with special care, they 
recover the formula obtained in the references |l^ from Feynman diagrams 
at lowest order in Zia. In our language, their identity reads: 



b,p,q 



i{q)TF{q,-p)v{p) 



ZlQ<l 



b,p,q 



'i-{i)TB{q,-p)v{p) 



ZlQ<l 



(55) 



However, one should refrain from trying to give a general interpretation to 

^*See also the section ^ of the present paper, which highUghts the differences between the 
retarded and Feynman propagators. 

The reason why scattering in non-relativistic quantum theory can be studied by solving 
the wave equation (even if the potential is non-static) is due to the fact that the Feynman 
and retarded free propagators are the same in a non-relativistic theory: 



•(P) = G°^"-{P): 



- p2/2m -I- ie ' 



(54) 



because the propagator associated to Schrodinger's equation has a single pole, at a positive 
energy 
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this result. Indeed, this identification works only at lowest order in Z\. However, 
the integrand in the Eq. (^l]) for ot is known to all orders in both Z\ and 
because the Dirac equation can be solved exactly with retarded boundary 
conditions. On the contrary, it is impossible to calculate all the corrections in 
Z\ in the perturbative calculation of ||, |l^ . We refer the reader to the sections 
H and 7.1 of the present paper, where we show that the retarded and time- 
ordered amplitudes are qualitatively different when considered to all orders in 
Z\ and Zi. In particular, the retarded amplitude has only a finite number of 
eikonal phases, while the time-ordered one can contain arbitrarily many of them. 
Therefore, it should be clear that the quantities calculated in ||l|, |[ |^ and in 
[H cannot be the same in general. 

The reason why Eq. (^5|) works at lowest order in Z\ is the following: if one 
considers the time-ordered amplitude only at lowest order in Z\, it simplifies 
dramatically into an expression that has at most two eikonal phases. If one 
focuses on the eikonal phases in u{q)Tp{q, —p)v{p) (see later the Eq. (106) for 
Tp, with p^,q^ > 0), one can show that (see p^, appendix E): 



u{q)TF{q,-p)v{p) 



-ieA{x^) _ 1] + ^p+^eA{v±) - l] + ^g-^eA{x±) _ ^ ^^+ieA{y±) -j^j 



(56) 



where the dots represent factors that are not needed for the argument. Doing 
the same thing for u{q)Tii{q, —p)v{p) (see Eq. (|75|), truncated at lowest order 
in Zi), one would get: 



J{q)TR{q, -p)v{p] 



ZlQ<l 



'ieA{xj 



-ieA{y± 



(57) 



where the dots represent exactly the same factors as in Eq. (^6|). One can see 
that the phase factors in the retarded amplitude differ from those in the time- 
ordered amplitude only by a global phase. It happens that this phase drops out 
when one is taking the modulus squared of those quantities and integrating over 
the momenta of the leptons, which explains Eq. (^5|). Because Eq. (^5|) relies 
on these properties of eikonal phases, it cannot be true for higher orders in Zi. 

Nevertheless, the calculation of reference is very interesting, because it 
indicates that the integration in Eq. ( ^ ) should be handled with great care. In 
particular, it seems that Coulomb corrections to ctt survive the integration over 
impact parameter b in the model defined by Eqs. (Ej) and (^. 

In view of the experimental result of |40| which observed an almost 

exact Z'^ scaling of the positron yield (or, equivalently, the number of produced 
pairs) apparently incompatible with large Coulomb corrections, one should men- 
tion a possible shortcoming of such a model. As explained in reference |12[, 



the Coulomb corrections in the integral Eq. (51) come entirely from the point 
fcj^ = 0. Therefore, they survive in gt only if there is no cutoff preventing zero 
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momentum transfers (either intrinsic to a theoretical calculation performed with 
a finite 7, or coming from the experimental setup). 

4 Propagators in the presence of one nucleus 

Even if the seeming discrepancy of the two approaches is now explained, there 
is still a paradox remaining. The work of ||, ||, Q indicates that ctt (or, 
equivalently, the retarded propagator) can be calculated exactly in the case of 
the collision of two ultra-relativistic nuclei. On the contrary, the perturbative 
expansion of |l^ leaves little hope that ui (or the Feynman propagator) can 
be calculated exactly. Why are these two propagators, calculated in the same 
background field, so different? This is the question we address in the next two 
sections, by showing that in the ultra-relativistic limit there is a simplification 
that allows to express in closed form the retarded propagator, but does not help 
to calculate the Feynman propagator. 

Let us first consider the case of the background field of only one nucleus 
of charge Z . It happens that for this case, the perturbative expansion of the 
scattering matrix can be summed into a closed expression when the nucleus is 
ultra-relativistic. This result will later appear as a building block in the expres- 
sion of the propagator in the case of two nuclei. For a nucleus moving at the 
speed of light in the positive z direction, the Fourier transform of the potential 
has the generic form A'^(fc) = 8{k~)v^k(h\}i. For a point-like Coulomb inter- 
action, A(fcj^) ^ Zefej^"^, but the following discussion does not depend on a 
specific form for A(A;^). Therefore, one could possibly take into account effects 
like an electromagnetic form factor for the nuclei. 

Instead of the propagator itself, it is simpler to deal with the scattering 
matrix T(q,p) (we reserve the calligraphic letter T for scattering matrices in 
the presence of two nuclei) defined by 

Giq,p) EE (27: fS{p - q)G"ip) + G°iq)T{q,p)G%p) . (58) 

The term of order n in the perturbative expansion of this object is given by 

T^{q,p)^{~ier j ^2^5(fcr)A(fci^)... j ^27r5(fc-)A(fe„^) 

xV+G"(p + fci>/+ • • • G\p + fci + • • • + K-i)tf+{2T:)H{p + fci + ■ ■ • + fc„ - g) . 

(59) 

At this stage, one can use the 5{k~) to perform for free all the integrations over 
the components. Using ^+ (j/ + r7i)?/+ = 2p^y+, we obtain 

T^{q,p) = [-ier2Ti5{p- - q-yj+ 
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X 5 -2 , (60) 

P'^ + ^t P+ + k+ + ■■■ + fc+-i - ^ + 

where is e > in the case of the retarded propagatoi]^, and sign(p~)e in the 
case of the Feynnian propagator, and where we denote: 



2 _ 2 



{p^ + ki^ + --- + kai_f . (61) 



The next step is to perform the integrals over the kf . A convenient trick is 
to introduce new variables Ai = + {ujf_-^ — ujf )/2p^ (and luq = 2p+p^) and 
write this integral as 



nl ^ J 2tt J 2tt ' 2p- 

perms, of the Ai 



Ai+iep H + 



(62) 



Note that the 1/n! is exactly compensated by the smii over permutations of 
the Ai (because these are interchangeable integration variables); it has been 
inserted for later convenience. Then, the following combinatoric formula 

^ i _ J i_ Ai H h A„ 



o-ee. 



where ©„ is the permutation group of [1, • • • , n], makes the various integrations 
almost independent (they are now coupled only by the 5 function). At this stage, 
one begins with the 5 function to get rid of A„ , and then performs successively 
the integrals over A„-i . . . Ai in the complex plane. The final answer for is 
extremely simple: 

/,! ^ ™^ ^ (64) 



n 



The transverse integral (second line of Eq. (|6^)) factorizes completely in the 
space of transverse coordinates, so that we obtain 

In^J d^x^ [A(cCi)]"e*(«^'P^)-"^ , (65) 

where A(a;_L) is the inverse Fourier transform of A(fc^). 
Collecting all the pieces, we have 

T„(,,p) . 2.Sip- - ,-V,sign(.,) / ,2^ J-»esign(.,)A(.,)]" ^,^^.^^,,^ 

(66) 



^For an advanced propagator, ep would be — e < 0. 
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and summing over n from 1 to +00 to get the fuU T{q,p) is now trivialp]: 



T{q,p) ^2tt6{p -q y+sign(ep) / d'^x 



-iesign(cp)A(a;^) _ -j^ 



(67) 

We can now be more specific, and write expUcitly the scattering matrix for 
the retarded (e^, > 0) and for the Feynman (sign(ep) = sign(p^)) propagators^ 



TR{q,p)^27rS{p- - q-)if+ d^x^ 



-?'eA(a:j_ ) 



^i{q±-p±)-x±_ 



(69) 



TF{q,p)^2nS{p -q >/+sign(p ) d^x± 



-ie sign(p )A(a:j_ ) 



J{q±-p±)-XA. 



(70) 

We observe that the retarded and Feynman results differ only by a sign(p ) 
appearing in two places. This sign will turn out to be essential when we go to 
the case of two nuclei, basically because the interaction with the second nuclei 
can change the sign of the p~ of the electrorP^ 




Figure 2: Schematic representation of the scattering of an e~ by one ultra- 
relativistic nucleus. The black dot denotes the scattering matrix T. The thin 
lines are free propagators, and the thick line is the trajectory of the nucleus. 



Using both Eq. ( [Sq ) and the above results for the scattering matrix T, the 
interpretation of the result is rather straightforward: the propagator is the 



'For the scattering matrix associated with a nucleus moving in the —z direction, replace 



P ,1 



,v+ by p+,q 



+ 



and use the appropriate A(a;^). 



^^For the advanced prescription, the result is: 
7A(<?,p) = -27r<5(p- -g-V+ jd^X 



(68) 



^^One can already notice here that it is the retarded version of the scattering which appears 
in the solution of the Dirac solution ^, ^, ^ This should not be a surprise, since 

solving the Dirac equation with initial boundary conditions involves naturally the retarded 
propagator, jll) noticed the sign{p~), but attributed to some mistake its absence in the 
solution of the Dirac equation. The above considerations show that the solution of the Dirac 
equation is correct, since it is Tji and not Tp that should appear in the solution with this 
type of boundary condition. 



20 



sum of two terms; one of them is the free propagator G° (nothing happens to 
the electron), and the second term contains the scattering matrix sandwiched 
between two free propagators. In space-time, the support of T is on the hyper- 
plane t = z, where the field of the nucleus is non-zero, and we can represent the 
term G^TG^ by the diagram of figure ^ 

The formulae of Eqs. (|69|) and dT^) also illustrate the discussion of section 



3.3, Indeed, they show that the retarded and Feynman prescription lead to the 
same result in the field of a single nucleus (i.e. in a background field that can 
be made static by a change of frame) if po > 0. This is perfectly consistent with 
the fact that the difference only shows up in a time-dependent background field. 



5 Propagators in the presence of two nuclei 
5.1 Formal derivation: Watson's series 

Now that we have an exact result in the case of one nucleus, we present a formal 
solution for the case of two nuclei that uses the previously obtained scattering 
matrices T{q,p) as its building blocks. This approach allows us to derive some 
results regarding the retarded and Feynman propagators in presence of two 
nuclei, that can be checked directly from perturbation theory. 

Let us assume that we have to solve some generic Lippmann-Schwinger equa- 
tion G = G'^ + G^VG (in this section, we use very compact notations for the sake 
of brevity; the previous formula is in fact an integral equation), and that the 
potential V receives contributions from a number of different scattering centers: 

a 

Let us assume also that the scattering matrices Ta for the individual scattering 
centers are known. They satisfy 

Then, the full propagator G resulting from the action of all the scattering 
centers can be formally written in terms of the T^. Indeed, it is a pure matter 
of algebra to check that the following object 

G = G° + ^G"r„G„ (71) 

a 

is a solution of the full Lippmann-Schwinger equation (see pages 750-752), 
provided the Gq satisfy 

G„ = G° + ^ G°r^G^ . (72) 

This formal solution amounts to a reorganization of the initial perturbative 
expansion, which rcsums infinite subsets of terms corresponding to the T^. De- 
spite this achievement, the problem of finding G in closed form is far from being 
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solved, because the equations ( |7^ ) are coupled integral equations that are very 
difficult to solve. The expansion in powers of the Tq, that emerges naturally 
from them is known in the literature as Watson's series. 

5.2 The case of moving scattering centers 

The set (^2|) of coupled equations has a very intuitive interpretation in the 
context of wave propagation. Indeed, the "Lippmann-Schwinger" equation ■0 = 
(j) + G^Vip (0 being the incoming wave, being a propagator for the free wave 
equation) is solved exactly in the same way by 

a 

V'a = + ^ G°r0V/3 ■ (73) 

In this solution, one can interpret G'^Taipa as the partial wave scattered by the 
scattering center a, and therefore is the wave seen by the scatterer a. The 
second equation then tells that the wave seen by the center a is made of the 
initial wave, and of contributions coming from the waves scattered by all the 
other centers. This interpretation is illustrated in figure |^ in the case of two 
centers. On this classical example, it is also intuitive that multiple reflections 




incoming 
wave 



Figure 3: Illustration of Watson's series in the case of two centers. The scattered 
waves bounce back and forth between the two scattering centers. 

of the wave on the two centers cannot happen if the centers are moving away 
at a speed larger than the velocity of wave propagation. 

However, the latter remark is correct only in the case of a non-relativistic 
wave equation. In the case of a relativistic wave equation, the limit where the 
two centers are moving at the speed of light is more intricate. In the next 
subsection, we show that in the case of the Dirac equation, a similar conclusion 
holds for the retarded propagator, but not for the Feynman propagator. 
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5.3 Retarded vs. Feynman propagators 
5.3.1 Retarded propagator 

If G" is the free retarded propagator, a typical contribution to the fuh retarded 
propagator is represented in the figure Q (left) when the two nuclei move at a 
speed V smaller than the velocity of light. Only this type of contribution can 
contribute to the retarded propagator, because the free retarded propagator 
connecting each scattering matrix is vanishing outside of the forward light-cone. 
This property has also important consequences in the case where the two nuclei 
are flying away at the speed of light, because the separation of two points lying 
respectively on the hyper-planes z — ±t is space-like, except if those points 
have z's of the same sign. This restriction forbids terms with more than two 
insertions of the scattering matrices (see figure ^, right). 




Figure 4: Typical contribution to the retarded propagator in the case of two 
nuclei. The thin lines are free retarded propagators, and the black dots are 
scattering matrices Tn. Left: the two nuclei are non-relativistic (velocity v < 1). 
Right: ultra-relativistic case; terms with more than two dots are forbidden 
(dashed line). 

In the perturbative expansion of the retarded propagator, this simplification 
of the ultra-relativistic limit appears in the following way: terms where inter- 
actions with the second nuclei are inserted between interactions with the first 
nuclei, like in 

TirG j^T2rG j^Tiji , (74) 

are vanishing {Tm and T2R are the retarded scattering matrices associated to 
the nuclei 1 and 2 respectively, in the retarded prescription (given by Eq. (|6^))). 
Indeed, it is immediate to check that such a term would have all its poles on the 
same side of the real-energy axis. More generally, the only terms allowed do not 
alternate interactions with the two nuclei: they can at most contain one "packet" 
of interactions with one nucleus, followed by a "packet" of interactions with 
the other nucleus. Using the above symbolic notations (integrations over the 
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momenta exchanged with the nuclei are implicit), the full retarded propagator 
in presence of two ultra-relativistic nuclei reads: 



i? — (-r_R 

+GrTibGrT2rG% + G%T2RG°jiTiiiG% . (75) 

These are the first three orders in the expansion of Eq. ( |7^ ) . All the following 
terms vanish because they involve factors like Eq. (|7^). This result is precisely 
the object that appeared in the solution of the Dirac equation |p], p|, |4| , 
confirming the fact that this approach in fact derived the retarded propagator. 
That was to be expected given the boundary conditions used to solve the Dirac 
equation. 

5.3.2 Feynman propagator 

Unlike the free retarded propagator, the free Feynman propagator can connect 
any pair of points in space-time. Therefore, the ultra-relativistic limit docs not 
forbid any contribution to the Feynman propagator. Physically, the additional 
terms correspond to the creation of additional e+e~ pairs, which are annihilated 
later so that only one electron is present when t —>■ +00. This is illustrated in 
figure H 

t z=t 





(2) \ 




/ Z 






z=-t 



Figure 5: Example of a contribution to the Feynman propagator that does not 
vanish in the ultra-relativistic limit. The thin lines are free Feynman propaga- 
tors, and the black dots are scattering matrices Tp. Point (1 ): a pair is created; 
point (2): a pair is annihilated. 

In perturbation theory, the main difference with the retarded case is that 
the analogue of Eq. (Q): 

TipGpT2FGpTip , (76) 

is not vanishing. This can be traced back into the fact that the free Feynman 
propagator has poles on both sides of the real-energy axis, and this in turn is 
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related to the fact that the free Feynman propagator can connect points with 
any time-ordering. As a consequence, the Watson's series for the Feynman 
propagator is infinite even in the ultra-relativistic hmit: 

/^O 

F — '-'F 

-\-L;rplipL;rp + pl2F'~^ F 

I /^O 'F /^O I /^O rp rp >^0 

-^^jr pllF'^ p^2F'-J F I Lrpl2F'~'F^lF'-^F 

I /^O /^O i^O I /^O /-^O rp >-^0 i^O 

-|-C_rj^i IFLt^ J 2F<-^F-' lf""^F + yj p-L2F'^ F^IF'^ f^'^F'^ F 

+ ■■■ . (77) 

One can note that the difference between the retarded and Feynman prop- 
agators is a feature specific to the relativistic nature of the Dirac equation. 
Indeed, for a non-relativistic wave equation like Schrodinger's equation, the 
retarded and Feynman ie prescriptions give the same propagator because the 
propagator has a unique pole, which has a positive energy (see footnote p^ . 
This remark is in agreement with the fact that the surviving terms involve the 
creation/annihilation of additional pairs, a purely relativistic effect. This inter- 
pretation in terms of pair creation/annihilation also highlights why the retarded 
and Feynman propagators are equally simple in the case of one nucleus, whereas 
they are not in the case of two nuclei: this is due to the fact that pairs cannot 
be produced by a single nucleus. 



6 Unitarity 

6.1 Calculation of | (Oout|0in) P 

So far, we have said nothing about the factor (Oout|0in) which appeared in the 
expression of the pair production amplitude, besides the fact that this factor is 
not just a phase. In this paragraph, we derive an expression for this vacuum- 
vacuum amplitude, that depends only on the Feynman propagator. In order to 
calculate | (Oout|0in) we start from the well known fact that the corresponding 
amplitude is the exponential of the sum of vacuum- vacuum diagrams : 

(OoutlOin) = e'^ , (78) 

with 




where the coefficients 1/2, 1/4, • • • are the symmetry factors of the the corre- 
sponding diagrams. Then, the probability of vacuum to vacuum transition is 
given by 

|(0out|0in)P =e^'^"^^ . (80) 
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To proceed, one could use cutting rules^ in order to compute the imaginary part 
of the vacuum-vacuum diagrams. However, this approach makes cumbersome 
the tracking of symmetry factors. For this reason, it is much simpler to remark 
that V is also given byp^ 

iV = Tt lii{l + ie4G°p) , (81) 

where the symbol Tr denotes a trace on Dirac's indices, as well as a trace over 
space-time. It is trivial to check that the expansion of the logarithm generates 
the series of diagrams in V, with the correct symmetry factors and a minus sign 
for the fermion loop. To begin with, one must write the Lippmann-Schwinger 
equation for the full Feynman scattering matrix^ (in presence of the two nuclei) . 
Formally, it reads 

Tf = -ie4 - ie4-GpTF = -ie4 - TFGpie4 , (82) 

from which we deduc^^ 

Tf + T* - -T*{G% + GV)Tf = ~T*{p^+^ + p^-'))Tf , (83) 

where we denote p^^\p) = 27r6'(±po)(2^+ m)S{p'^ — m^). 

Noticing then that closed loops of retarded propagators are zero, and that 
G% =G%,~ p^~\ we have 

= Tr \n{l + ie4^G%) 

= Tr ln(l -I- ie^G^ - ie^p^''^) 

= Tr ln(l + ie4G"F + (1 + ieiG°F)TF p^^'^) 

= Tr ln(l + Tpp'^^^) . (84) 

From there, one obtains 

-2lmV = i{V -V*) =-Tr ln((l + Tj.>(-))(1 + Tf/?'"))) 

= -Tr in((i + rj?p(-)rV"' + (Ti. + r;)p(-))) 

= -Tr ln(l - T^p^+^FP^'^) ■ (85) 



Therefore, we have 



^"in order to use cutting rules, it is essential to notice that, since the sign of the energy 
flowing in the photon lines is not fixed, one cannot exclude contributions where the cut divides 
the diagram in more than two connected pieces. This is precisely what happens in terms like 
the second diagram of figure H 

^^This formula, together with Eq. wq), is well known for the vacuum-to- vacuurn amplitude. 
In the literature on pmi production in nuclear collisions, it already appears in |42j (see their 
Eq. (45)). However, [M does not work out its modulus squared | (Oout|0in) P- 

22 Let us recall that Gf = G% + G^TpG^. 
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From this formula, we know (Oout|0in), up to an irrelevant phase. This expression 
shows that the vacuum-vacuum amplitude is completely determined once we 
have set up some approximation scheme that gives the Feynman propagator. 
The very reason for this is the nature of the Lagrangian of Eq. (^. Indeed, this 
Lagrangian implies that Wick's theorem holdsQ so that we can express the full 
2n-point Green's functions in terms of Gp only. 



6.2 Probability to produce n pairs 

Since there is some confusion in the literature regarding unitarity in ultra- 



relativistic heavy ion collisions ll 14 13, 16, 43, M ttSL W6 , it is useful to 



derive the probability of production of n pairs. The corresponding transition 
amplitude is 

Mn{{p^,q^}) = {e+{pl)■■■e+{pn)e-{ql)■■■e-{qn)out\Oin) . (87) 

This amplitude can be related by a reduction formula similar to Eq. (^) to 
the Fourier transform of the amputated 2n-point time-ordered correlator (with 
appropriate spinors for the final states). Again, we pull out the factor (Oout|0in), 
in order to get rid of vacuum-vacuum diagrams in the second factor: 

/n in \ (Oout|TV;(?/i)---^(yn)V'(xi)---Vj(a;„)|Oin) _ 

(OoutlOin) ^ r-T = (0out|0i„)G2„({x,,,?/J) . 

\Uout|Uin/ 

(88) 

At this stage, it is important to notice that since the Lagrangian Eq. (g) does 
not contain any dynamical field that couples to the fermions. Wick's theorem 
applies to the second factor G2ni{xi, yi}) so that we can write it in terms of the 
full Feynman propagator: 

G2n{{xi,yi}) ^ ^ e{a)GF{xi,ya{i))- ■ ■GF{Xn,ya(n)) , (89) 

where e{a) is the signature of the permutation a (required when permuting 
fermion fields). This property extends to the amputated correlators and their 
Fourier transform, so that we can write directly the transition amplitude asp^ 

M„({p,,qJ) = (OoutlOin) fM["(9'T(i))XF(g^(i),-pi)w(pi)]--- 

• • • mqa{n))TF{qain), -Pn)v{p„)] , (90) 

where Tp is the Feynman scattering matrix in presence of the two nuclei. 

^■'This would not be true if we had kept the photon kinetic term F^^F^'-'^ / i in the La- 
grangian. 

^^This formula is equivalent to the Eq. (11) of_jk2j. However, these authors do not make 
use of it to calculate the exact Pn (see our Eq. (pj|)). 
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The integrated probability of producing exactly n pairs is then (we have 
changed pi —^ —pi in the second line) 



(OoutlOi, 



E <<')i^^T\[H(l^)'^F{q^,P^W^ + nl)T^{p,,q,^e,)u{q,(^^)] (91) 



where p^^\p) = 2TT9{±pn)6{p^ - m^). We can see now that we obtain closed 
chains like: 

tr [Ti^ (gi,Pi)(2/i +m)Tp{pi,q2)(g(2 + (92,^2) (2^2 + ™)'^f (P2, gi)('/i + m)] . 

(92) 

For instance, for n — 2, cycles come in two sizes, illustrated on figure |[ A 
systematic tool to construct these loops is the (unique) decomposition of per- 
mutations in products of disjoint circular permutations. It is a trivial matter 
of combinatorics to find that the number of p-cycles in (S„ is n\/p{n — p)!, and 
that the number of permutations that are made of ai 1-cycles, 02 2-cycles,. . .,a„ 
n-cycles (ai + 2a2 + • • • + rta„ = n) is: 



n! 



ai!---a„! I'^i 



(93) 



The signature of such a permutation is e{a) = (— 1)" HiC"!)"'- With the fol- 
lowing compact notations for "links"rj and "loops" : 



L{p',p) 



TrL" = 



d'^q 
(2^ 

d^P. 



n 



\ (27r)4 

the total probability to produce n pairs is 



V{p',q)p^+\q)TF{q,p)p'^-Hp) 



tr[L{pi,p2)L{p2,P3) ■ • •L(p„,pi) 



(OoutiOi„)p j2 n 

ai +2a2H \-nan —n i—1 



(-1)"' /TrL'' 



v * 



(94) 



(95) 



We should emphasize here the fact that the individual P„ arc functions of the 
Feynman scattering matrix Tp, while the average number of pairs n = nP„ 
has a simple expression in terms of the retarded scattering matrix Tr (see 
Eg. (ig|)). 

^^Note that L = in a background field that cannot produce pairs, like the field of a single 
nucleus for instance. 
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6.3 Another look at unitarity 

We can now check unitarity without making any approximation, by using the 
previous formula for the P„, and Eq. (|8^) for the vacuum- vacuum transition 
probability: 



n— n— ai 4-2a2 H \-nan—n i—1 



Pn =|(0out|0i„)P^ n 



^ g-Tr ln(l-r;p( + )r^.p<-))gTr In(l-L) ^ ^ ^gg^ 

The compensation of the two exponentials is obvious from the first of Eqs. (Q). 

It therefore clarifies the long standing "unitarity problem" floating around in 
the literature related to pair production in nuclear collisions |l^, |l^, |l^ 
^ The problem can be stated as follows: perturbation theory seems 

to lead to cross sections that are too large in order to comply with unitarity (or 
to production probabilities larger than 1). It appears that the factor (Oout|0in) 
(the modulus of which is smaller than one) has been overlookec^ in the litera- 
ture: only the connected piece of Feynman diagrams has been considered, but 
not the disconnected vacuum-vacuum diagrams. In this paper, we have shown 
that the factor (Oout|0in) naturally emerges from the reduction formula for pair 
production (and therefore that calculating only the connected diagrams leads 
to an incomplete result), and that this factor restores unitarity. 

Moreover, since the factor (Oout|0in) depends only on the Feynman prop- 
agator, the above considerations provide a way to make approximations that 
preserve unitarity. It is sufficient to use the same approximate Gi? in the calcu- 
lation of the Pn and in the calculation of the vacuum-vacuum amplitude. 

For such an approximation scheme to be consistent, we have also to verify 
that it gives positive P„. Indeed, the fact that the sum of Eq. ( p6| ) is 1 leaves 
open the possibility that some probabilities could come out negative. To prove 
that all the probabilities are positive (and hence smaller than one, because 
their sum is one), it is convenient to introduce a "generating function" for the 
probabihties P„: 

F{x) = e-T'- '"(1-^)6^^ ^"(1-^^) , (97) 

such that -F(l) = 1 and P„ = P'^"-'(0)/n!. Then, one notices that it can be 
rewritten as 

F{x) = e-T'- ln(l+«')eTr ln(l+xttt) ^ ^gg^ 

with t{q,p) = u{q)Tp{q, p)v{p) . The operator tt^ is positive, and we can write 
the generating function in terms of its eigenvalues , which are positive: 

^(-)-n^- (99) 

i 

^^Two exceptions are Ji2| and [[^st, which noticed that unitarity is related to the vacuum to 
vacuum amplitude. However, the authors of |45| calculated this amplitude by requiring that 
unitarity is preserved (in addition, they only did that in an approximation that leads to a 
Poissonian distribution). One needs an independent jailculation of (Oout|0in) in order to claim 
that this factor restores unitarity. The authors of [Ua] give an exact expression of {Oout|0in) 
equivalent to our Eq. (|8l|), but unfortunately do not exploit it to check unitarity. 
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We therefore see that all the coefficients in the Taylor expansion of F{x) are 
positive independently of the approximation made for Tp, which proves that all 
the Pn are positive in this framework. 

6.4 Nature of the multiplicity distribution 

One can also note that Eq. (pi) for the probability of producing n pairs is not a 



Poisson distribution. This fact contradicts the solution proposed in |16, ^ for 
the unitarity problem, where a Poissonian distribution is obtained. However, 
one can see from Eq. (|9^ ) the nature of the approximation that would lead to 
such a distribution: a Poisson distribution is obtained if one drops all the TrL" 
for n > 1. Indeed, this drastic (and a priori not justified) simplification leads to 

P„-e-™^^^, (100) 
n! 

a Poisson distribution of average — TrL. This explains the somewhat confusing 
statement of jl^, ^ Q saying that the "perturbation theory prediction" 
for the probability to produce one pair (i.e. Pi in which one would forget the 
factor I (Oin|0out) that is —TrL in our notations) should not be interpreted as a 
probability (because it can be larger than 1) but instead should be reinterpreted 
as the average number of pairs produced in a nuclear collision. Modulo the 



approximation of Eq. (100), we agree with this statement, except for one thing: 
—TtL is not the perturbation theory prediction for the probability to produce 
one pair; perturbation theory applied correctly indicates that this probability is 
— I (Oout|0in) P Tri < 1. In fact, the present analysis makes the following clear: 
even if Tp (and hence L) were known exactly, using —TrL for the probability to 
produce one pair would still violate unitarity. The unitarity problem does not 
really come from perturbation theory, but from forgetting the contribution of 
the vacuum-vacuum amplitude. 

Thanks to the generating function introduced above, it is easy to find what 
the average number of produced pairs isp^: 

iP„ = F'(l) = -Tr[L(l - L)-i] Tri" ■ (1^3) 



+ 00 +00 

n=0 



By the same method, we can find the variance of the number of pairs: 

^-n^ = -Tr[L(l - L)-^] - Tr[L2(l - L)-^] . (104) 

^*One can chock that this formula is equivalent to Eq. (po|), by using the Lippmann- 
Schwinger equations for Tp and Tr. One can obtain successively: 

Tf = [1 - TRp^-^r^Ta , 

L = [i- T^p(-)rM?3p'+>w-'][i - w-)]-^ , 

[1 - Tr>(-'][i - W-)] - [r^>(+)Tflp(-)] = 1 , (101) 

and finally 

n = -Tr [L{1 - L)-^] = -Tr [T^p'+'W-'] , (102) 
which can be cast into Eq. ( |5o| ) thanks to(/+m = Jospin andj/— m = X]spin ^(P)^(P)- 
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Since n ^ n? — v? ^ the exact probability distribution cannot be a Poisson 







(Triy 






Tr (L^) 



Figure 6: Two types of pairings in the probability to produce two pairs. Even 
if the amplitude factorizes because we neglect photons connecting directly the 
various fermion lines, the cross-section does not. 



distribution. This is a consequence of the permutations in the final state when 
squaring the amplitude, which can lead to assign to the same pair^the electron 
of one fermion line and the positron from another fermion line. This is illustrated 
in figure |^ for the case of two pairs, where one can see clearly the origin of the 
terms (TrL)^ and TrL^. One also sees that the term in Ti correlates the 
emission of the two pairs, and therefore prevents the probability distribution 
from being Poissonianp"]. This explains the fact that a Poisson distribution is 
obtained only if one neglects the TrL" for n > 1. 

The authors of |Q attempted to size the departure from a Poisson distri- 
bution in the Magnus model, and found it to be around 1%, which may justify 
for most practical purposes to calculate only n (which can be calculated ex- 
actly according to our analysis) and plug it in a Poisson formula for P„. This 
is the approach followed in started from the lowest order (only 

one photon is exchanged with each nucleus, and (Oout|0in) is approximated by 
1) Pi in the background field of two nuclei, "reinterpreted" as n. followed 
the same approach by starting from the formula given in Eq. (53) of ^J, which 
was assumed to be the exact formula for the pair production amplitude. If this 
assumption were correct, its square should be the exact probability to produce 
one pair, and should therefore be smaller than 1. This is not the case as noted 
in who "reinterpreted" this quantity as the average number of pairs pro- 
duced in a collision in order to save unitarity. These authors did not realize that 
their remark was directly pointing to the result we have proven in our section 
^: that the retarded solution of the Dirac equation does not give the pair pro- 
duction probability, but rather the average number of pairs (without any need 
to "reinterpret" anything). 

^^The problem comes from the fact that particles are produced in pairs here. As an example 
of model in which particles are produced individually, one can consider a toy model in which 
a scalar field (f>{x) is coupled to a background classical source j{x) by Ci^t = i(^)</'(^)- This 
model is exactly solvable, and the production probability is found to bePoissonian (see j23|, 
pages 163-170). It is also possible to understand from this remark why [ha] obtained a Poisson 
distribution. Indeed, this paper modeled e+e~ pairs as elementary fields (quasi-bosons), so 
that the problem of pairings in the final state never shows up. 

'^''Another argument against a Poisson distribution is that it would contradict the fact that 
n can be calculated exactly, while the individual Pn cannot. 
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7 Strategies of approximation for Gp{x,y) 



7.1 Limit Zia < 1 

Since the series in Eq. (|7^) giving the fuU Feynman propagator cannot be 
summed exactly, one must use approximations in order to simphfy it. We 
present in this section two different approximation schemes which lead to closed 
form expressions for the Feynman propagator. 

The simplest approximations one can think of are truncations of the Watson 
series. Such a simplification is obtained in a natural way if one assumes that 
only one nucleus has a large electric charge. If Zi is small enough so that we 
have 

Zia < 1 - Z2a , (105) 

then we can neglect most of the corrections in Zia. Since we are interested in 
pair production, the kinematics requires at least one interaction with the first 
nucleus (otherwise, the pair cannot be produced on-shell). We are therefore go- 
ing to keep only terms with at most one insertion of . In that approximation, 
Tip is just — ie^i, and we drop any term with two or more insertions of Tip in 
the Watson series. It is immediate to verify that the only terms left in Gf are 

p ~ Lt^ 

+G°,(-i6^i)G°, + G%T2fG''p 
+G%{-i(4i)G%T2fG% + G%T2pG"pi-^di)G% 
+G%T2pG''pi-ie4i)G°pT2pG°p . (106) 

This approximation for Gp is in fact the starting point used by ^, ^ in their 
approach to the problem of pair production. The terms that participate to pair 
production (i.e. having at least one interaction with each nucleus) are displayed 
on figure 1^. The term labeled (1) in this figure corresponds to the creation of an 
on-shell electron and an off-shell positron. The positron is subsequently put on 
its mass-shell by an additional interaction. For the process (2), both the electron 
and the positron are first created off-shell, and then interact independently to 
go on-shell. 

This approximation, which leads to an analytically tractable expression, suf- 
fers however from several limitations. The obvious one is that for the collision 
of two nuclei like gold for RHIC, the charge of the nuclei is not small enough to 
justify this approximation, and the inequality ( |l05| ) is not satisfied^. Moreover, 
this approximation leaves out some terms that seem physically important. In 
particular, it does not include any term where both the and the e~ in the 
pair interact with the two nuclei. Such a term would indeed have four scattering 
matrices. 

^^For the case Zi ~ Z2, a term has to be added were the roles of Zi and Z2 are exchanged. 
The accuracy of the result obtained via this procedure is claimed to be better than the percent 
level for RHIC energies ||, 
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Figure 7: Dominant pair production mechanisms in the approximation Z^a <C 
1. The thin hnes are free Feynman propagators, the black dots are scattering 
matrices, and the open dots denote interactions at lowest order. 



7.2 Limitation in the number of intermediate pairs 

Another way to approximate the Feynman propagator is to start from the defi- 
nition _ 

GF{x,y) = jT — |— ^ . (107) 

Then, one can use a completeP^ set of states {|riin)}ri=o - oo counts the number 
of pairs, other continuous indices hke the momenta of the particles have not been 
written explicitly) and insert in the previous equation an identical operator 
constructed as 

Tl = 

Separating the contribution from ri = 0, we obtain 
GFix,y)= (Oin|T^(y)V(a:)|Oin) 



E'^°° (Ooutj^in) 
(OoutIO 



+ L n n'" (nin|TV;(y)^(x)|Oin) ■ (109) 



In this formula, n can be seen as the number of extra pairs produced (and 
then destroyed) in the course of the evolution of the electron. Another way to 
simplify the Feynman propagator is therefore to truncate the previous sum, and 
keep only the term obtained with n = 0, that is 

GF(x,y)«(Oi„|T^(y)V(x)|Oin) • (110) 



•^^ Strictly speaking, these states form a basis of the subspace containing only Fock states 
with zero electric charge. This is sufficient here since they will be contracted with the vacuum. 
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One can note that both the exact Feynman propagator Gp and the approx- 
imation (0in|T7/;(?/)^/'(a;)|0in) are solutions of the Dirac equation {i0x — eJ^{x) — 
m)G{x,y) = d{x — y). Their difference is therefore a solution of the homoge- 
neous Dirac equation. It is in fact immediate to verify that all the terms with 



n > 1 in the right hand side of Eq. (109) are solutions of the homogeneous Dirac 
equation. 

The calculation of section [3.2.2| in fact gives the answer for this correlator in 
terms of retarded and advanced propagators only: we have an exact expression 
for the Fourier transform G++(g,p) of (Oin|Ti/)(j/)V'(a;)|Oin), which reads 



G++{q,p) - GR{q,p) + GAiq,p) + Gs{q,p) 



(111) 



with Gs{q,p) given in Eq. (^). We see that contrary to the exact GF{q,p), 
this correlator can be written in closed form, in terms of the exactly known Gr 
and Ga'- 



G++{q,p) = - Gniq^p) + GA{q,p) 



1 
2 

(2^ 



2^5{e - m^) G%iq)TRiq, k){H + m)r^(fc,p)G°i(p) , (112) 



where Tr and Ta are the retarded and advanced scattering matrices in presence 
of the two nuclei. Each of them can contain up to two of the 1-nucleus scattering 
matrices Tr^a (see Eq. ([ts])). 



Z=t 




z=-t 



Figure 8: JVew contribution to pair production in the approximation Gp ~ 
(Oin|Ti/)(y)V'(a;)|Oin)- The thin Unas are free propagators, the black dots are 
scattering matrices. 



A common feature of the two expressions Eqs. (B^) and (111) and Eq. ( |112[ ) 
for the approximation of Gp by <^Oin|T^/;(2/)'0(a;)|Oinyis that they contain terms 
involving up to four scattering matrices on individual nuclei. Therefore, this 
approximation contains a little more than the previous one (based on Zia ^ 1), 
since it includes a selected subset of the terms of order four in the Watson series. 
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The physical meaning of those terms is the following: an off-shell pair is first 
created. Then, the electron and the positron independently scatter off each of 
the two nuclei. No additional pair is created or annihilated. These terms are 
illustrated in figure |[ 

8 Conclusions 

In this paper, we have studied various aspects of the problem of pair produc- 
tion in the collision of ultra-relativistic heavy ions, focusing in resolving the 
discrepancies of the existent literature. 

By showing that the inclusive cross-section of pair production is related to 
retarded amplitudes, while the Feynman amplitude gives the exclusive cross- 
section of single pair production, we found that the discrepancy between two 
methods used to attack the problem of pair production lies at a deeper level 
than expected: these two methods do not calculate the same physical quantity, 
and their result should not be compared directly. 

Then, we have studied the propagator of an electron in the electromagnetic 
background field created by two colliding nuclei. It appeared that the exact 
retarded propagator can be obtained in closed form if the nuclei are ultra- 
relativistic. However, such a simplification does not occur for the Feynman 
propagator, which can only be expressed as an infinite series. This problem 
arises in any background field that can produce pairs of particles, and can 
be traced back into the relativistic nature of the wave equation governing the 
electron field. This observation indicates that in the collision of two ultra- 
relativistic nuclei, the inclusive cross-section of pair production can be obtained 
exactly, but not the more exclusive ones. Experimentally, it would therefore be 
desirable to measure both types of cross-sections. 

Within the model of |^, the inclusive cross-section, which is expressed in 
terms of the exactly known retarded amplitude, seems to contain Coulomb cor- 
rections as well. However, the calculation of completely integrated cross-sections 
is questionable in this model. Indeed, the problem raised by Lee and Milstein 
occurs at zero momentum transfer (or at infinite impact parameter), precisely 
where the model is expected to break down. In particular, any infrared cutoff 
on the momentum transfer, coming either from an improved theoretical model 
or from experimental cuts, could drastically reduce these corrections. 

We have also analyzed the unitarity puzzle, and shown that pair production 
probabilities satisfy all requirements of unitarity if the factor | (Oout|0in) P is 
correctly taken into account. In addition, since everything depends on the 2- 
point Feynman Green's function, unitarity is preserved if one starts from an 
approximation of this propagator. A side product of this analysis is that the 
multiplicity distribution is not Poissonian. 

Finally, in addition to the approximation Zia ^ 1 ^ Z2 used in [To| , 
we have presented a completely different approximation for Gpix^y) that leads 
also to closed expressions, and seems to contain more of the relevant physics for 
pair production by two heavy nuclei. 
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